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Abstract. We discuss the relevance of the scalar modes appearing in chiral theories with sponta- 
neous symmetry breaking such as the NJL model for nuclear matter studies. We show that it depends 
on the relative role of chiral symmetry breaking and confinement in the nucleon mass origin. It is 
only in the case of a mixed origin that nuclear matter can be stable and reach saturation. We de- 
scribe models of nucleon structure where this balance is achieved. We show how chiral constarints 
and confinement modify the QCD sum rules for the mass evolution in nuclear matter. 



Introduction 

One crucial question of present day nuclear physics is the interplay between the 
nuclear many-body problem and the nucleon structure where both chiral symmetry 
breaking and confinement are involved. In this talk we will develop a viewpoint where 
the nuclear attraction is associated with the in-medium fluctuation of a (chiral invariant) 
background scalar field which is itself at the origin of the constituent quark mass. 
Nuclear stability is ensured with the incorporation of the nucleon response to this scalar 
field. This response depends on the quark confinement mechanism inside the nucleon. As 
we will see this question can be rephrased in terms of the relative weight of spontaneous 
chiral symmetry breaking and confinement in the origin of the nucleon mass. 



Chiral nuclear matter description including the effect of confinement 

As a starting point let us consider the relativistic mean-field approaches initiated by 
Walecka [1] where the nucleons move in an attractive scalar and a repulsive vector 
background fields. It provides an economical saturation mechanism and a well known 
success is the correct magnitude of the spin-orbit potential where the large vector and 
scalar fields contribute in an additive way. Now the question of the very nature of these 
background fields and their relationship with the QCD condensates has to be elucidated. 
To address this question we formulate an effective theory parametrized in terms of 
the fields associated with the fluctuations of the chiral condensate in a matrix form 
= O + it ■ %) by going from cartesian to polar coordinates, i.e., going from a linear 

to a non linear representation : ^ = a + ix ■ % = SU = (f K + s) exp yvt ■ <t> n / ' f-^j . In 

ref. [2] we made the physical assumption to identify the chiral invariant scalar field, s = 



S — fx, associated with radial (in order to respect chiral constraints) fluctuations of the 
condensate, with the background attractive scalar field. The normalization is such that 
the symmetry breaking piece of the Hamiltonian is given by H x $b = — J d^rF n M\ o. It 
follows that the evolution of the condensate is given by : 

{qq) vac F n 2F* F^M* 

The first equality is valid to leading order in the pion field and the second equality holds 
for a dilute system. The quantity s = (s) is the expectation value of the nuclear scalar 
field. In this picture the nuclear medium can be seen as a shifted vacuum characterized 
by a chiral order parameter S = F n + s and the nucleon mass will depend in some way 
on the scalar field M N = M N (S). The dynamics of the scalar field is encoded in the 
chiral effective potential, W(s) = (M^) e ffS 2 /2 + associated with chiral symmetry 
breaking. Notice that the effective sigma mass has no reason to coincide with the mass 
of the /o(600) [3, 4]. In the low density limit the expectation value of the scalar field 
is s = —gap I iMcs)eff where g a = (dM^ / ds) s=0 is the scalar coupling of the nucleon. 
The pion-nucleon sigma sigma commutator receives a contribution from the scalar field 

and one from the pion cloud, On = o Na ■ 

4T iGn) = F*go ojf iond0 ^ = [ d\ \m\ (N\ 4> 2 |A0 . (2) 

The quark condensate evolution can be compared with the one of the nucleon mass 
which, for a dilute system, is : 

^ x , J_ { _ So^ti Q Na P l _ ga_ Q Na P (3) 
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We see that this result deviates from the Ioffe sum rule [5] which predicts that the quark 
condensate and the nucleon mass have the same evolution with density controlled by the 
pion nucleon sigma term [6, 7]. Firstly the pion cloud piece of the sigma commutator, 
^{pioncloud) ^ contr [butes to the condensate evolution. It does not contribute to the mass 
evolution, otherwise chiral constraints would be violated (such as the presence of a term 
in m K in the NN potential, forbidden [8] by chiral symmetry). In fact its influence on 
the mass vanishes in the chiral limit and hence it is a small effect which we ignore. 
Only the scalar piece, o N ° pion , should then enter the mass evolution. There is however 
another deviation from the Ioffe sum rule which is the presence of the weighting factor, 
geF K /MN, where the scalar couping constant, g a , depends on the nucleon structure. 
Only in case where the nucleon mass entirely originates from chiral symmetry, such as 
in the linear sigma model, or naive additive NJL model (g CT = M N /F K ), the Ioffe sum 
rule is recovered but corrected from pionic effect. The pion-nucleon sigma term entering 
these evolutions is an important piece of experimental information. It is obtained from 
the Feynman-Hellman theorem : ajy = m{dMx/dm) ~ 50MeV. According to previous 
works [8, 9, 10], we expect for the pion cloud contribution O ^ wncloud ^ ~ 20MeV. 
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The numerical value ot the non pionic piece has to be a^ a =0^ — 0^ ~ 
50 — 20 ~ 30MeV. We will come back to this important constraint later on. 

The Hartree energy density of nuclear matter (including omega exchange) writes in 
terms of the order parameter s as : Eq/V = £o = / (4d 3 p/ (2n) 3 ) &(pf — p) E*(s) + 

W(s) + g 2 v /2m 2 p where E*(s) = \J P 2 + Mff(s) is the energy of a nucleon with 
momentum p. Here two serious problems appear when the nucleon mass has a pure 
symmetry breaking origin : M^(s) = Mn + g a s. The first one is the fact that the chiral 
effective potential, W(s), contains an attractive tadpole diagram which generates an 
attractive three-body force destroying matter stability [11]. The second one is related 
to the nucleon substructure. According to the lattice data analysis of Thomas et al 
[9], the nucleon mass can be expanded according to Mjv(m|) = ciq + ci2tn 2 l + ^fflj + 
L n (m n , A) + where the pionic self-energy is explicitly separated out. While the ci2 
parameter is related to the non pionic piece of the kN sigma term, a$ is related to the 
nucleon QCD scalar susceptibility. The important point is that {a^)i att ~ — 0.5GeV 3 
[9] is essentially compatible with zero in the sense that it is much smaller than in 
our chiral effective model, {®A)chiral = 

-3F„g a /2(M*) ef f ~ -3.5 GeV 3 , where the 
nucleon is seen as a juxtaposition of three constituent quarks getting their mass from 
the chiral condensate [12]. The common origin of these two failures can be attributed 
to the absence of confinement [13]. In reality the composite nucleon responds to the 
nuclear environment, i.e., by readjusting its confined quark structure [14]. The resulting 
polarization of the nucleon is accounted for by the phenomenological introduction of the 
positive scalar nucleon response, K N $, in the nucleon mass evolution : M N (s) = M N + 

goS + \ Ktvss 2 + This constitutes the only change in the expression of the energy 

density but this has numerous consequences. In particular the parameter is modified : 
CI4 = {ciA)chiral i} ~ f ■ The value of C = (F 2 /2M N ) k n $ which reproduces the lattice 
data is C~ 1.25 implying a strong cancellation effect in a 4 [12, 10]. Moreover the scalar 
response of the nucleon induces an new piece in the lagrangian Jz? s 2 NN = — K^ss 2 NN/2 
which generates a repulsive three-body force able to restore saturation [15]. 

The restoration of saturation properties has been confirmed at the Hartree level [12] 
with a value of the dimensionless scalar response parameter, C, close to the value 
estimated from the lattice data, taking g a = M^/F K and g m adjusted near the VDM 
value. The next step has been to include pion loops on top of the Hartree mean- field 
calculation [10, 16]. One possibility is to use in-medium chiral perturbation theory 
but we prefer to use a standard many-body (RPA) approach which includes the effect 
of short-range correlations (g r parameters fixed by spin-isospin phenomenology) and 
A — h excitations. The calculation which also incorporates rho exchange has no real free 
parameters apart for a fine tuning of C (around the lattice estimate) and of g m (around 
the VDM value). The inclusion of spin-isospin loops improves the quality of the result, 
in particular they bring down the compressibility close to the accepted value. We stress 
the relatively modest value of the correlation energy (~ — lOMeV), much smaller than 
what is obtained from iterated pion exchange (planar diagramm) in in-medium chiral 
perturbation theory. This effect is mainly due to the strong screening of pion exchange 
by short-range correlations. For completeness we mention that 

we have also performed a full relativistic Hartree-Fock calculation with the notable 
inclusion of the rho meson exchange which allows to also reproduce asymmetry proper- 



ties of nuclear matter. Details can be found in ref. [13]. 



Status and relevance of the nuclear scalar field 

Although the concept of a scalar field has been widely used for nuclear matter studies 
[1] its precise origin or meaning is still a controversial subject. The problem is that there 
is no sharp scalar resonance which would lead to a simple scalar particle exchange. In 
our approach instead we stress the chiral aspect of the problem. As soon as we start 
from a model which gives a correct description of chiral symmetry breaking in the QCD 
vacuum such as the Nambu-Jona-Lasinio model (NJL), the emergence of a scalar field 
linked to the quark condensate cannot be avoided. For this purpose we introduce the NJL 
model in the light quark sector limited to describe the main mesons : the pion, the sigma, 
the rho, the a\ and the omega mesons. The lagrangian is : Jzf = \[f (iy^d^ — m) \y + 



(Gi/2) [(v/v) 2 + (V//75T v/) 2 j - (G2/2) [(xfrft yY + {ytY5?vY + {wtw) 
It is generally accepted that this model gives an excellent description of vacuum chiral 
symmetry breaking, with the appearence of a constituent quark mass and a soft Gold- 
stone pion mode. In addition in the so called delocalized version the inclusion of a form 
factor, depending on a cutoff, A, at each quark leg of the interaction generates a mo- 
mentum dependence of the quark mass in agreement with lattice calculation. However 
the approach is notoriously not fully satisfactory due to the lack of confinement : in 
particular unphysical decay channels of vector mesons in qq pairs may appear. In that 
respect it has been shown that adding a confining interaction on top of the NJL model 
solves the problem of unphysical qq decay channels of mesons with masses larger than 
twice the constituent quark mass [3]. In addition the scalar meson arising around twice 
the constituent quark mass is pushed at higher energy by confinement and the c(600) 
comes out as a broad %% resonance. It has also been demonstrated that the effect of the 
confining interaction, although crucial for the on-shell properties of vector and scalar 
mesons, has only little influence for the low momentum physics relevant for the nu- 
clear many-body problem [4]. Hence, our attitude will be to derive an effective low 
momentum theory where the mass parameters for scalar et vector mesons will not be 
the on-shell mesons masses but simply mass parameters associated with the inverse of 
the corresponding correlators taken at zero momentum. Technically this can be done by 
rewriting the NJL lagrangian in a semi-bozonized form and integrating out quarks in the 
Dirac sea using a path integral formalism. The physical meaning is simply a projection 
of qq vacuum fluctuations onto meson degrees of freedom. Keeping only the relevant 
terms for nuclear physics purpose, the resulting low momentum effective lagrangian has 
the form : 

J2f = -f^d^S c d^ c -W(y = g0S S) + -F 2 Ml — tr f (U+Ui-2) 

-^{co^co^ + v^-v^ + a^-a^). (4) 



5? is a chiral invariant scalar field whose vacuum expectation value coincides with 
the constituent quark mass M and S c is the canonical one. We also introduce an 
effective scalar field (5%// = (F n /Mo)y = F n + s normalized to F n in the vacuum 
which coincides with the chiral invariant S field introduced in the second section and W 
is the corresponding chiral effective potential. The matrix U is U = exp(<&/F) where 4> 
is the orthoradial pion field and the parameter F can be identified with the pion decay 
constant parameter F n . The quantity & = (Mq/F 71 )(F 71 + s) is the expectation value of 
the scalar field which is expected to decrease in the nuclear medium. The various Fs 
functions are standard NJL loop integrals which also enter the mass parameters My 
and Ma for vector and axial-vector mesons. These masses also depend on the vector 
coupling constant G 2 . In particular the ratio g v /M v (gy being the quark- vector coupling 
constant) is fixed to G 2 . Finally the pion mass emerges as M\ = mMo/GiF 2 . We have 
a priori four parameters, G\ , G 2 , the cutoff A and the bare quark mass m. We use : A = 
1 GeV, m = 3.5MeV, G\ = 7.8 GeV 2 . We thus obtain for the vacuum quark mass 
at zero momentum : Mq = 31\MeV and (qq) = —(2S6MeV) 3 . We actually constrain 
G 2 to be close to the VDM value (G 2 ) VDM = g v /M v = (2.65/0.770) 2 GeV~ 2 . For 
the nuclear matter calculation we take G 2 = 0.78 (G 2 ) VDM F n = 93.6MeV, M K = 
137.8MeV. With this set of parameters the low momentum mass parameters are : 
(M a )eff = 653MeV, M v = \256MeV andM A = 1398MeV. 

In the pure NJL picture, at finite baryonic density, the value of the constituent quark 
mass, which is the expectation value of the scalar field 5?, is modified. It can be obtained 
self-consistently from a gap equation modified by the presence of a Fermi sea. However 
in the real world baryonic matter is not made of independant constituent quarks but of 
clustered objects, the nucleons. These nucleons are embedded in the scalar background 
field, and the nuclear medium can be seen, as said before, as a shifted vacuum. The 
nucleon mass will depend in some way on the scalar background field and the energy 
density of symmetric nuclear matter at the Hartree level reads : 



Eo 
V 

(5) 

The expectation value for the scalar field is self-consistently obtained by minimization 
of the energy density, 

which constitutes an in-medium modified gap equation with (M^(S^)) being the 
nucleonic scalar density. The scalar coupling constant of the nucleon to the effective 
scalar field (which is normalized to F K in the vacuum) is : 



( \ <&\ M ofdM N \ 
(Ss)eff(y) = 1 r[-^) 



which depends crucially on nucleon structure. For instance if the nucleon mass fully 
originates from confinement (bag models), dM N /dy = 0, the scalar field just decouples 




FIGURE 1. Left panel: Mass of the quark (dashed line), of the diquark (dot-dashed line) and of the 
nucleon (full line) versus the relative deviation of the scalar field with respect to its vacuum value. Right 
panel: Binding energy of nuclear matter versus nuclear matter density in units of normal density. The full 
line corresponds to the full result and the dashed line represents the Hartree result. The dot-dashed line 
corresponds to the contribution of the Fock term and the dotted line represents the correlation energy. All 
the numerical inputs are given in the text. 



from the nucleon, ((gs)eff = 0)- In this case there is no shift of the vacuum and the scalar 
field is thus an irrelevant concept for nuclear matter studies. On the other extreme if the 
nucleon mass fully originates from chiral symmetry breaking (naive additive NJL, chiral 
soliton), then the nucleon mass in the medium is affected by the scalar field associated 
with the dropping of the chiral condensate. However in this case the attractive tadpole 
destroys stability. Only in the case where the nucleon mass has a mixed origin, the 
scalar background field can contribute to the nuclear attraction without destroying the 
stability and saturation properties. In that case by rearranging its quark structure linked 
to the confinement mechanism, the nucleon reacts against the scalar field generating 
effectively repulsive three-body forces. The origin of this repulsion lies in the decrease 
of the scalar coupling constant of the nucleon. In short a possibly important part of 
the saturation mechanism is associated with the progressive decoupling of the nucleon 
from the scalar field associated with the dropping of the chiral condensate. In the next 
section we will introduce nucleon models capable of achieving the balance between 
large enough attraction and sufficient reaction. Of course, one falls here in the modelling 
uncertainties. However there is that a stringent constraint for the numerical value of the 
scalar nucleon coupling constant which is model dependent, from the value of the free 
nucleon sigma commutator. 



Effect of confinement: simple models for the in-medium nucleon 

We now come to the last point of this paper, namely the modelling of the nucleon 
mass origin and the scalar response of the nucleon defined from the second derivative 
of the nucleon mass with respect to the scalar field. For a nucleon made of the simple 
adjonction of three NJL constituent quarks (or a NJL quark and a NJL diquark) the scalar 
coupling constant is independant of the scalar field and there is no scalar response. The 



importance of the response is related to the respective roles of chiral symmetry breaking 
and confinement in the generation of the nucleon mass. In a previous work [15] we 
have introduced a model of a nucleon made of three constituent quarks bound together 
by a confining harmonic force. The magnitude of the scalar response which followed 
was too small to prevent the collapse of nuclear matter. We will come back later to 
this type of model. A possibility of improvement is to reduce the relative role of chiral 
symmetry breaking. This can be achieved by considering a nucleon made of a quark and 
a sufficiently light diquark to leave enough room for confinement. A practical advantage 
is that a three-body problem is transformed into a simpler two-body problem. Beside 
this simplification, there are theoretical and phenomenological reasons to favor a quark- 
diquark model of the nucleon with relatively light scalar-isoscalar diquark. For instance 
the work of Shuryak et al on hadronic current-correlation functions based on a random 
instanton vacuum [17] finds a strong attraction in the scalar-isoscalar channel leading to 
a diquark with a mass about 400 MeV . 

As discussed in a set of works of Bentz et al (see ref. [11]), it is possible to construct 
from the NJL model a nucleon which has a diquark component. The diquark mass is 
also medium dependent since it depends on the constituent quark mass. If we take 
for the the coupling constant in the diquark channel, G\ = 0.92 Gi, we obtain for the 
vacuum diquark Mo = 398.5 MeV which turns out to be nearly equal to the constituent 
quark mass. In [11], a nucleon scalar response is obtained through the inclusion of 
an infrared cutoff Hr ~ 200MeV in the Schwinger proper time regularization scheme. 
Such a prescription implies that quarks cannot propagate at relative distance larger than 
1/flR, hence mimicking a confinement mechanism. Here we propose to incorporate 
confinement in a more direct way. Since the diquark is in an anti-triplet color state, 
it is physically plausible that a string develops between the quark and the diquark as in a 
QQ meson. We thus introduce a confining potential between the quark and the diquark: 
V(r) =Kr 2 /2. In the non relativistic limit, the problem reduces to solve the Schrodinger 
equation for a particle with reduced mass /i, placed in an harmonic potential. In this limit 
the mass of the (in-medium) nucleon is given by : 



M Ni y ) =M ( y ) + M Di y ) + \^ ) with „ = »*L (7) 

Taking for the string tension a standard value K = (290MeV) 3 , we obtain for the 
vacuum nucleon mass Mm = \30AMeV . The nucleon mass origin splits roughly into 
a chiral symmetry breaking component (60%) and a confinement component (40%). 
The vacuum value scalar coupling constant of this nucleon to the effective scalar field 
is g a = (gs) e ff(Mo) = 7.14. This leads to the value of the non pionic piece of the 

sigma term calculated in the model : o^° pion) = F K g a (Ml/ (Ml) eff ) = 30MeV, as 
was required. In the linear sigma model where g a = Mn/F k = 10 we recover the Ioffe 
sum rule just corrected from pionic effects. With confinement the value of the scalar 
coupling constant is reduced with respect to the linear sigma model or additive NJL 
model (ga = Mn/F k = 10) and the mass evolution is slower than the condensate one 
(see eq. (1,3)). The suppression of the pionic contribution to the mass evolution further 
accentuates the difference between the mass and condensate evolutions. For instance at 
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FIGURE 2. effective scalar coupling constant versus the relative deviation, <t> = (Mo — M) /Mq = \s\ /F K , 
of the scalar field with respect to its vacuum value for the linear confining potential (dashed line), the 
quadratic linear potential (dotted line) and for the quark-diquark model (full line). 



normal nuclear density the condensate has dropped by ~ 30%. With the value g a ~ 7 
deduced above, the mass reduction is significantly lower, ~ 13%. 

In order to show that such a model is capable of describing the saturation properties 
of nuclear matter we calculate the energy of symmetric nuclear matter in the Hartree 
approximation, using eq.(5, 6). The resulting curve displays a saturation mechanism 
driven by the scalar nucleon response (Kjvs, proportional to the second derivative of 
the nucleon mass with respect to ^) which has a positive value. Said differently the 
scalar coupling constant, dM^/d^, is a decreasing function of \s\ or the density. This 
translates into the fact that the nucleon mass stabilizes or even increases with increasing 
|*| (see fig. 1). However the binding is nevertheless not sufficient unless we decrease 
artificially the vector coupling constant G2 at a value much smaller than the VDM result. 
In order to improve the description, although this is not necessarily consistent with our 
present nucleon model, we add on top of the Hartree mean field result the pion loop 
(Fock term and correlation energy) contribution obtained in our previous work [10]. 
Taking the value of G2 at the value quoted previously, G2 = 0.78 (G2) VDM , we obtain a 
decent saturation curve shown in fig. 1. Likely a fully consistent calculation within the 
model of the pion loop energy would modify the result but a fine tuning on G2 would be 
presumably sufficient to recover the correct saturation curve. The lesson of this simple 
model calculation seems to confirm our previous conclusions. The confinement effect 
(scalar response of the nucleon) is able to stabilize nuclear matter and the pion loop 
correlation energy helps to get the correct binding energy. 

We have shown that an acceptable quark-diquark model of the nucleon makes plau- 
sible the role of the background scakar field in the nuclear binding. It is interesting to 
investigate if other confining mechanisms can achieve the same result. For this we have 
also studied models where the nucleon is made of three constituent quarks moving in 
a mean-field linear confining potential but shifted with a constant attractive potential 
mimicking short range attraction : V = [(1 +yo)/2] (K2r — 2Vq). This model has been 
successfully utilized for baryon spectroscopies studies by Jena et al [18]. We do not aim 
to justify this particular equally mixed scalar and vector confining potentials, the main 
motivation being the existence of analytical solutions. The energy of the lowest orbit, 



solution of the Dirac equation, is 



E(M) =M-2V + ^[K 2 x q with^ solution of x 4 q + 2 ° xt - (2.3381 1) 3 = 

and the mass of the in-medium nucleon (in absence of CM correction) is M^{^ = M) = 
3E(M). Hence the quark mass contribution (essentially the chiral symmetry breaking 
contribution) to the quark orbital energy and then to the nucleon is reduced due to the 
presence of the attractive shift, —2Vq, leaving more room for the confining part. It is 
possible to show that the scalar coupling constant (still omitting CM correction) can be 
written as (gs)eff(^) = 3 (Mq/F k ) q s where q s = J d 3 r (u 2 — v 2 ) (r) is the quark scalar 
charge. We see that the scalar field contribution to the sigma term is represented by 
the usual integrated scalar quark density as in bag models. In practice we also include 
in the numerical calculation the effect of CM correction using the results quoted in 
ref. [18]. If we take K 2 = (300MeV) 2 and V = 200MeV it is possible to obtain a 
saturation curve but the saturation has the tendency to come too early. Certainly this 
point deserves a more detailed study. Here we wish to concentrate on the main result, 
namely a decreasing scalar coupling constant when increasing \s\ as demonstrated by the 
dashed curve on fig. 2. We also checked that replacing the linear potential by a quadratic 
potential, V = [(l + y )/2] (K 3 r 2 - 2V ) with 2K 3 = (300MeV) 3 , V = 200MeV, one 
obtains similar results as depicted on fig. 2 (dotted curve). It is worthwile to notice that 
this model differs from the one used in [15] by the introduction of the constant attractive 
shift —2Vq. Again this shift allows to reinforce the role of confinement in the origin of 
the nucleon mass. Also shown on fig. 2 is the behavior of the scalar coupling constant for 
the quark-diquark model. In this case, the decrease at low density is less strong which 
translates into a softer equation of state. According to a preliminary study based on 
a variationnal relativistic calculation the strong dropping beyond \s\ /F K « 0.2 (which 
roughly corresponds to normal density) might be to some extent an artefact of the non 
relativistic approximation. 
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